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Abstract 

FRW universe in Horava-Lifshitz (HL) gravity model filled with a combination of dark mat¬ 
ter and dark energy in the form of variable modified Chaplygin gas (VMCG) is considered. The 
permitted values of the VMCG parameters are determined by the recent astrophysical and cos¬ 
mological observational data. Here we present the Hubble parameter in terms of the observable 
parameters P^mO) ^vmcgOi Hq, redshift z and other parameters like a, A, 7 and n. From Stern 
data set (12 points), we have obtained the bounds of the arbitrary parameters by minimizing the 
test. The best-fit values of the parameters are obtained by 66%, 90% and 99% confidence levels. 
Next due to joint analysis with BAO and CMB observations, we have also obtained the bounds 
of the parameters {A, 7) by fixing some other parameters a and n. The best fit value of distance 
modulus is obtained for the VMCG model in HL gravity, and it is concluded that our model 
is perfectly consistent with the union2 sample data. 

Pacs no : 04.60.Pp, 98.80.Qc, 98.80.-k 


1 Introduction 

For several decades we have tried to give birth to a well constructed quantum gravity theory that can reconcile 
with the general theory of Relativity. With the discovery of the late cosmic acceleration [n 1 a n [5] at the 
turn of the last century, the unified theory became all the more necessary. Several footsteps can be found in 
literature that aims at producing a UV complete theory. Motivated by the success of the theory proposed 
by Lifshitz [5] in solid state physics, Horava proposed a gravity theory widely known the Horava-lifshitz (HL) 
gravity [7]. Taking the UV limit into account, HL gravity has a Lifshitz like anisotropic scaling as t —>■ Pt and 
x® —>■ Zx® between space and time. As this is characterized by the dynamical critical exponent z = 3, it breaks 
the Lorentz invariance, while in the infra red limit, the scale reduces to z = 1, i.e., it is reduced to classical 
general relativistic theory of gravity in the low energy limit. Even if we add a ^ term, ‘a’ being the scale factor, 
with the Friedmann equation [SI Hi El m min HU we will have HL gravity equations. 

In the original proposal, Horava 015] considered projectability condition with/without the detailed balance. 
The detailed balance condition was originally proposed in order to reduce the number of operations in the action 
(i.e., number of independent coupling constants) and to simplify some properties of the quantum system. As a 
result the form of potential in the 4-dimensional Lorentzian action is restricted to a specific form in terms of a 
3D Euclidean theory. But from cosmological point of view, this condition leads to major obstacles and hence 
should be abandoned. 

On the other hand, the fundamental symmetry of the theory namely the foliation-preserving diffeomorphism 
invariance leads to the projectability condition. In particular, one should have 3D spatial diffeomorphism and 
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the space-independent time reparametrization. As the lapse function is essentially the Gauge degree of freedom 
associated with the time reparametrization, so it should be independent of space coordinates. This is termed as 
projectibility condition. Due to this projectibility condition, the Hamiltonian constraint is not a local equation 
satisfied at each spatial point but an equation integrable over a whole space, i.e., a global Hamiltonian constraint. 
Note that absence of projectability condition and imposing a local version of the Hamiltonian constraint lead 
to phenomenological obstacles and theoretical inconsistencies. 

However, the important drawback of projectibility condition is that : due to the global Hamiltonian con¬ 
straint the general relativity could not be recovered from the Horava gravity for arbitrary small lambda. Also 
the projectibility condition from the point of view of condensed matter physics may not be appropriate for 
describing the (quantum) gravity. 

In extragalactic astronomy, type la supernovae have a characteristic light curve. The similarity in absolute 
luminosity profiles of nearly all known type la supernovae helps astrophysicists to treat them as a standard 
candle. Using this fact the observations Riess et al and Perlmutter et al in 1998 mm led us to the conclusion 
that our universe is going through a cosmic acceleration. There are two popular ways that provide theoretical 
support to the present day accelerated expansion of our universe. One is to modify the geometric part of 
Einstein equation, leading to the ideas of modified gravity and the other is to consider the universe to be filled up 
uniformly by some exotic matter possessing negative pressure. The Friedmann equation ^ = —AttG (p -I- 3p) /3 
requires the condition {p + 3p) < 0 for accelerated expansion (d > 0). As density is an ever positive physical 
quantity, we see the Equation of state (EoS) parameter must be negative and also less than —1/3. Such a 
negative pressure creating substance is aptly termed as dark energy (DE here after). Recent observational data 
indicates that DE occupies 73% of the whole matter-energy of our universe. Theoretically we can find many 
proposed DE candidates. Variable modified Chaplygin gas (VMCG) is one among them. The equation of state 
of VMCG is given as [El [TT] 


p = ap — 



( 1 ) 


Here we consider /3(a) = /3oa where /3o and 7 are constants. 

In 1904, Chaplygin [E] introduced Chaplygin gas (CG), whose EoS is given by 


p = - 


P 

P 


( 2 ) 


where /3 is a positive constant. An attractive feature of the model is that it behaves as dust-like matter at an 
early stage and as cosmological constant at later stages i.e. the CG behaves like a pressure less fluid for small 
values of scale factor and as a cosmological constant for large values of scale factors which tend to accelerate 
the expansion. At the beginning of this century, the CG model went through a series of modifications, all in 
the quest of a suitable candidate for DE. 

In |2Q], Bento et al for the first time gave the idea of Generalised CG having the EoS 


p = - 


P 


(3) 


Whereas in 2010, Jamil |2T] presented a model in which the new generalized Chaplygin gas interacts with 
matter. Benaoum in 2002 further modified CG model and proposed the Modified CG (MCG) model obeying 
the equation of state [32 HSl |34l [25] 


p = ap ■ 


P 


(4) 


with Of > 0 and 0 < n < 1. When a = |, this EoS shows radiation era at one extreme (when the scale 
factor is vanishingly small) while ACDM model at the other extreme (when the scale factor is infinitely large). 
At all stages it shows a mixture. Amidst these there also exist one stage when the pressure vanishes and 
the matter content is equivalent to pure dust. On further modification VMCG came into existence. In [18] . 
Jamil investigated the evolution of a Schwarzschild black hole in the standard model of cosmology using the 
phantom-like modified variable Chaplygin gas and the viscous generalized Chaplygin gas. 

Present day trend of literature says that the combination of DE with modified gravity together gives more 
interesting results. In 2010 Park showed in [E] that the Friedmann equation in Horava gravity contains 
additional a“^, and cosmological constant terms as the effective DE and he predicted that these terms may 
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be responsible for cosmic acceleration. These terms, coined as Horava effective DE, were tried to be constrained 
on a few occasions [271128] . The best ever approach was via phenomenological parametrization of the relevant 
physical quantities that have been well-studied. Park himself considered the widely used CPL parametrization. 
B. C. Paul and his coleagues have studied constraints for different exotic fluid model in the back ground of HL 
gravity [211 HO] • In this work, we shall study the limits of the DE parameters constrained by different data 
sets. First we will set the constraint for the closed universe and analyze it. Later on we deal with the range of 
parameters for open and flat universe. It is quite understood that any observational constraints on HL gravity 
do not enlighten the discussion about the well-known conceptual problems and instabilities of the theory, nor 
it can address the questions concerning its validity. Therefore in the present analysis the HL gravity has to be 
considered as a phenomenological model. The same holds for the Chaplygin Gas model as well. 

The paper is organized as follows: The basic generalized equations for HL gravity are given in section (2). 
Various dimensionless density parameters have been discussed in section (3). The main mechanisms which will 
be followed to analyze the data is briefly given in section (4). Lastly, a brief summary and a fruitful conclusion 
have been drawn in section (5). 


2 Basic Equations in Horava-Lifshitz Gravity 

It is convenient to use the Arnowitt-Deser-Misner decomposition of the metric which is given by dUlIllIS] 

ds2 ^ + gij {dx^ + N^dt) {dx^ + N^dt) . (5) 

Here, N is the lapse function, Ni is the shift vector, gij is the spatial metric. The scaling transformation of the 
coordinates reads as : t —>■ and a;* —>■ lx*. The HL gravity action has two constituents, namely, the kinetic 
and the potential term as 

Sg = Sk + Sy = J dtd^Xy/gN [L^ + L„) (6) 

where, the kinetic term is given by 


Sk = J dtd^Xy/gN 

where, the extrinsic curvature is given as 


2 - XK"^) 


_ 9ij 

“ 2N 


(7) 


( 8 ) 


The number of invariants, while working with the Lagrangian, can be reduced due to its symmetric property 
. This symmetry actually is known as detailed balance. Considering this detailed balance the expanded form 
of the action becomes 


Sg= dtd^Xy/gN 


2 - \K^) 


n'^g'^RijR^^ 




8(3A- 1) 


2a;4 

(1 - 4A)i?2 


2 u; 2 ^ 
+ AR- 3A^ 


(9) 


here is the Cotton tensor and all the covariant derivatives are determined with respect 

to the spatial metric, is a totally antisymmetric unit tensor, A is a dimensionless constant and k, oj and 

g are constants. 

Now to set the matter-tensor the total energy density and pressure, pm and Pm respectively will be taken 
which will contain in itself the impact of the baryonic density (pb), dark matter density (pdm), etc. 

Assuming only temporal dependency of the lapse function (i.e., N = N{t)), Horava obtained a gravitational 
action. Using FRW metric with iV = 1 , gij = a^{t)^ij , iV* = 0 and 

dr'^ 

"iijdx^'dx^ = -- + r^dVt\ 

1 — Kr^ 
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where K = —1, 1,0 represent open, closed and flat universe respectively and taking variation of N and gij we 
obtain the Friedmann equations [Ml [30] 


= 


6(3A- 1) 


(Pm + Pr) + 


6(3A- 1) 




+ 




1(3A- l)a4 8(3A- 1) 


8(3A- l)^a2 


( 10 ) 


H + 


3if2 


4(3A- 1) 


“b Pr'^r) 


4(3A- 1) 


8(3A- l)a4 


3kV"A2 ■ 
8(3A- 1)_ 


K^gi'^kK 
8(3A- l)^a2 


( 11 ) 


The term proportional to ^ is an unique contribution of HL gravity, which can be treated as “Dark radiation 
term” munis] and the constant term is the cosmological constant. 

The conservation equation of matter is 


Pm “b (Pm ~b Pm) — 0 

and that of radiation is 


( 12 ) 


Pr + 3iA (pr +Pr) — 0 (13) 

where, Gcosmo = i 67 r( 3 A-i) ■^ilh the conditions = 1 and Ggrav = 

Solving the DE conservation equation for variable modified Chaplygin Gas, the expression for DE density 
for VMCG will be 


Pvmcg — Pvmcgoi^ + [A + (1 — A)(l + z)^] “+^ . 

1 

' “I" 1 

, A = —, where C is the integration constant. 

PrmcgO 


where Pvmcgo = 


3(rt+l)Po 

7 


c 


3 Observational Constraints on EOS Parameters 


(14) 


Using Eqs. (lO)-(ll), the Friedmann’s equation can be rewritten as: 

SttG, . f A\ K 

= 


g '{Pvmcg + Pdm + Pr) + l 2Aa'^ 


• / 1 

H -\- ^ — AtiCt \ Vvmcg 

We define the following dimensionless density parameters: 
(i) for matter component: 


-( 


3A 

V4A^ ~ T 


K 


(ii) for curvature: 


(iii) for cosmological constant: 


o 


n - ^ 


On — 


772 


A 


Another dimensionless parameter for expansion rate is defined as: 

H{z) 


E{z) = 


Ho 


Using the above parameters, the Fridmann equation can be rewritten as: 

02 

E‘^{z) = ^vmcgoF{z) + O£;mo(l + z)^ + Oro(l + (1 + z)^ + Og + Oifo(l + z)'^ 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 
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where 


( 22 ) 


F{z) = Pvmcgo{l + [^1 + (1 - A)(l + Z^] ° + i 

Here is for radiation and fldmo is for dark matter. At the present epoch E(z = 0) = which leads to 


n 


vmcgO 


+ ^dmO + ^rO + 


02 

“iCO 

Hn 


+ Hq + ^KO — 1 


(23) 


4 Observational Data Analysis In HL Universe 

In this section, we perform a detailed observational data analysis [32131131131133 using Stern data. We also 
study the model under Stern+BAO and Stern+BAO+CMB joint analysis. The mechanism that we will use in 
the present work is the minimum test from theoretical Hubble parameter with the observed data set and 
find the best fit values of unknown parameters for different confidence levels (66%, 90%, 99%). In the table 
given below, we present the 3 column Stern data. 


z 

H{z) 

cr(z) 

0 

73 

± 8 

0.1 

69 

± 12 

0.17 

83 

± 8 

0.27 

77 

± 14 

0.4 

95 

± 17.4 

0.48 

90 

± 60 

0.88 

97 

± 40.4 

0.9 

117 

± 23 

1.3 

168 

± 17.4 

1.43 

177 

± 18.2 

1.53 

140 

± 14 

1.75 

202 

± 40.4 


Table 1: The Hubble parameter H{z) and the standard error a{z) for different values of redshift z. 


4.1 Stern {H{z)-z) Data Set 

In this sub-section our theoretical model of VMCG in HL gravity is analyzed, using the observed values of 
Hubble parameter at different redshifts (twelve data points) listed in observed Hubble data by Stern et al [36) . 
The observed values of Hubble parameter H{z) and the standard error a{z) for different values of redshift z 
are given in Table 1. A statistical hypothesis is proposed and its validity is tested at different confidence levels. 
For this purpose we first form the statistics as a sum of standard normal distribution as follows: 


\Stern 


E 


{H{z) - Hobs{z)Y 
0 -2(2;) 


(24) 


L = J e (25) 

where H{z) and Hobs{z) are theoretical and observational values of Hubble parameter at different redshifts 
respectively. Here, Hobs is a nuisance parameter and can be safely marginalized. Hq is the present value of 
Hubble parameter and its value is fixed at Hq = 72 ± 8 Kms“^ Mpc“^. From the DE model we see that the 
two most important parameters are A and B{a). Here we shall determine best fit value of the parameters {A vs 
7) by minimizing the above distribution x'stern fixing the other unknown parameters with the help of Stern 
data. According to our analysis the best fit values of A vs 7 is presented in Table 2. We now plot the graph 
for different confidence levels. Our best fit analysis with Stern observational data support the theoretical range 
of the parameters. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) contours are plotted in 
fig.l for a = 0.001, n = 0.1, = 0.01, = 0.02, ^ko = 0.01, Oq = 0.01, Vldmo = -03. 
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H(z)-z data (Stern) 
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Fig.l 


Fig.3 


Fig. 2 

The contours are drawn for 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) confidence 
levels. 

Fig.l shows the variations of A against 7 for a = 0.001, n = 0.1,njo = 0.01,12^0 = 0.02,0^0 = 
0.01, flo = 0.01, = -03 with respectively for the H{z)-z joint analysis. 

Fig.2 shows the variations of A against 7 for a = 0.001, n = 0.1,fljo = 0.01, Oj-o = 0.02, = 

0.01, flo = O.OljfldmO = .03 with respectively for the H{z)-z+BAO joint analysis. 

Fig.3 shows the variations of A against 7 for a = 0.001, n = 0.1,fljo = 0.01,flj.o = 0.02,0x0 = 
0.01, Oq = 0.01, = -03 with respectively for the i?(2:)-2;+BAO+CMB joint analysis. 


4.2 Stern + BAO Data Sets 


Here we resort to a joint analysis in the sense that the BAO peaks are incorporated in the stern data. The 
Baryon Acoustic Oscillation (BAO) peak parameter value was proposed by [37] and we shall use their approach. 
The pioneer as far as the detection of BAO signal is concerned, is considered to be the Sloan Digital Sky Survey 
(SDSS) survey. The survey directly detected the BAO signals at a scale ^ 100 MPc. The analysis that is 
followed is actually the combination of angular diameter distance and Hubble parameter at that redshift. This 
analysis is independent of the measurement of Hq and does not contain any particular dark energy. Here we 
examine the parameters A vs 7 for VMCG model from the measurements of the BAO peak for low redshift 
(with range 0 < z < 0.35) using standard analysis. The error corresponds to the standard deviation, where 
the distribution considered is Gaussian. We know that the Low-redshift distance measurements are very lightly 
dependent on different cosmological parameters, the EoS of dark energy and have the ability to measure the 
Hubble constant Hq directly. The BAO peak parameter is defined by 


A = 


F(2i)1/3 




(26) 


Here E{z) = H{z)/Hq is the normalized Hubble parameter, the redshift Zi = 0.35 is the typical redshift of the 
SDSS sample and the integration term is the dimensionless comoving distance for the redshift Zi. The value of 
the parameter A for the flat model of the universe is given by A = 0.469 ± 0.017 using SDSS data [33 from 
luminous red galaxies survey. Now the function for the BAO measurement is given as. 


Xbao — 


{A - 0.469)2 
(0.017)2 


The total joint data analysis (Stern+BAO) for the function may be defined by 


(27) 


Z _ Z I z 

\total XStern ' \BA0 


( 28 ) 
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According to our analysis the best fit values of A vs 7 for the joint scheme is presented in Table 2. Finally we 
generate the closed contours of A vs 7 for the 66 % (solid, blue), 90% (dashed, red) and 99% (dashed, black) con¬ 
fidence limits depicted in fig.2 for a = 0.001, n = 0.1, Oio = 0.01, Qro = 0.02, ^Iko = 0.01, Qq = 0.01, fldmo = -03. 


4.3 Stern + BAO + CMB Data Sets 


The angular scale of the first acoustic peak is measured through angular scale of the sound horizon at the 
surface of last scattering. This is one of the most interesting geometrical probe of dark energy. The information 
is encoded in the CMB (Cosmic Microwave Background) power spectrum. The definition of the CMB shift 
parameter is given in [31IM1I10]- It is not sensitive with respect to perturbations but are suitable to constrain 
model parameter. This is the property that we will use in this analysis. The CMB power spectrum hrst peak 
is the shift parameter which is given by 



(29) 


where Z 2 is the value of redshift at the last scattering surface. From WMAP7 data of the work of Komatsu et al 
m the value of the parameter was obtained as 77. = 1.726 ± 0.018 at the redshift z = 1091.3. The function 
for the CMB measurement can be written as 


2 (77- 1.726)2 

XCMB - ( 0 . 018)2 


(30) 


Now when we consider three cosmological tests together, the total joint data analysis (Stern-|-BAO-|-CMB) 
for the function is defined by 


XtOTAL — Xstern + XbAO + XCMB ( 31 ) 

Now the best fit values of (A, 7 ) for joint analysis of BAO and CMB with Stern observational data support the 
theoretical range of the parameters are given in Table 2. The 66 % (solid, blue), 90% (dashed, red) and 99% 
(dashed, black) contours are plotted in fig.3 for a = 0.001, n = 0.1, Oio = 0 . 01 , 0^0 = 0.02, = 0.01, Oq = 

0.01,Od^o = .03. 


Data 

A 

7 

Xmm 

Stern 

0.0512871 

-2.02163 

7.17375 

Stern + BAO 

0.0536315 

-2.03587 

768.149 

Stern + BAO + CMB 

0.942108 

13.3671 

9554.2 


Table 2: The best fit values of A, 7 and the minimum values of x^- 

4.4 Redshift-Magnitude Observations from Supernovae Type la 

The main evidence for the existence of dark energy was provided by Supernova Type la experiments. The 
existence of dark energy is directly related to the redshift of the universe. Therefore, since 1995 two teams of 
High redshift Supernova Search and the Supernova Cosmology Project have been working extensively, and in 
their effort they have discovered several type la supernovae at the high redshifts dJllllIll]. The observations 
directly measure the distance modulus of a Supernovae and its redshift z HUS]. Here we will consider the recent 
observational data, including SNe la which consists of 557 data points and belongs to the Union2 sample |44) . 
From the observations, the luminosity distance c?l(z) will determine the dark energy density which is defined 

by 


dhiz) = (1 + z)Ho f 

Jo 


dz' 


The distance modulus (distance between absolute and apparent luminosity of a distance object) for Super¬ 
novae is given by 
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/i(z) = 5 logio 


dL{z)/Ho 
1 MPc 


+ 25 


(33) 


The best fit of distance modulus as a function fi{z) of redshift z for our theoretical model and the Supernova 
Type la Union2 sample are drawn in fig.4 for our best fit values of the parameters. From the curve, we see that 
the theoretical VMCG model in HL gravity is in agreement with the union2 sample data. 


5 Discussions 

In this work, we have considered the FRW universe in HL gravity model filled with a combination of dark matter 
and dark energy in the form of variable modified Chaplygin gas (VMCG). We present the Hubble parameter 
in terms of the observable parameters ilio, ^ko, ^dmOj ^vmcgo and Hq with the redshift z and the other 
parameters like H, n, 7 and a. For these parameters we have chosen the specific numerical values consistent 
with observations. From Stern data set (12 points), we have obtained the bounds of the arbitrary parameters 
by minimizing the test. Next due to joint analysis of BAO and CMB observations, we have also obtained the 
best fit values and the bounds of the parameters (A, 7). We have plotted the statistical confidence contour of 
(A, 7) for different confidence levels i.e., 66%(dotted, blue), 90%(dashed, red) and 99%(dashed, black) confidence 
levels by fixing observable parameters VLdmO, ^r0^ ^ko, ^0 and Hq and some other parameters like, n and a, 
etc. for Stern, Stern+BAO and Stern+BAO+CMB data analysis. 

From the Stern data,the best-fit values and bounds of the parameters (A, 7) are obtained. The output values 
are shown in Table 1 and the figure 1 shows statistical confidence contour for 66%, 90% and 99% confidence 
levels. Next due to joint analysis with Stern + BAO data, we have also obtained the best-fit values and bounds 
of the parameters (A, 7). The results are displayed in the second row of Table 1 and in figure 2 we have plotted 
the statistical confidence contour for 66%, 90% and 99% confidence levels. After that, due to joint analysis with 
Stern+BAO+CMB data, the best-fit values and bounds of the parameters (A, 7) are found and the results are 
shown in Table 1. The figure 3 shows statistical confidence contour for 66%, 90% and 99% confidence levels. 
For each case, we compare the model parameters through their values and by the statistical contours. From the 
comparative study, one can get an idea about the convergence of theoretical values of the parameters with their 
values obtained from the observational data set and how it is altered for different chosen set of other parametric 
values. 

Finally the distance modulus /i(z) against redshift z has been drawn in figure 4 for the theoretical model 
of the VMCG in HL gravity for the best fit values of the parameters and the observed SNe la Union2 data 
sample. So the observational data sets are perfectly consistent with our predicted theoretical VMGG model in 
HL gravity. 

The present study discover the constraint of allowed composition of matter-energy by constraining the range 
of the values of the parameters for a physically viable VMCG in HL gravity model. In a nut-shell, the conclusion 
of this discussion suggests that even though the quantum aspect of gravity have small effect on the observational 
constraint, but the cosmological observation can put upper bounds on the magnitude of the correction coming 
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from quantum gravity that may be closer to the theoretical expectation than what one would expect. 
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